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In this paper we consider multimode multiphoton Jaynes-Cummings model, which consists 

O " 

' of a two-level atom, initially prepared in an excited atomic state, interacting with N modes 

bJ)' of electromagnetic field prepared in general pure quantum states. For this system we show 

that under certain conditions the evolution of the Wigner function at the phase space origin 
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provides direct information on the corresponding atomic inversion. This relation is also valid 
even if the system includes Kerr-like nonlinearity, Stark shift effect, different types of the 
initial atomic state as well as moving atom. Furthermore, based on this fact we discuss for 
the single-mode case the possibility of detecting the atomic inversion by means of techniques 



<**] ■ similar to those used for Wigner function. 

q-( PACS numbers: 42.50.Dv,42.50.-p 

> 

CO . I. INTRODUCTION 

IT) 

^ : n 

Jaynes-Cummings model (JCM) [lj has continued to be a subject of not only theoretical studies 
[ — . \ but also experimental investigation (, e.g. see [2|]). This model in the simplest form is described 

as a two-level-atomic system interacting with an electromagnetic field (for a review see, e.g. ^]). 



o 
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Many of quantum features of the JCM have been predicted and observed: among the most well 

5h ' 

known is the revival-collapse phenomenon (RCP) of the atomic inversion (a z (T)) 



a a. 



RCP 



is arised from the presence of multiple exchange of photons between the radiating atom and the 
cavity mode. Observation of RCP for (a z (T)) has been performed using the one-atom mazer Q], 
which is more sophisticated than the dynamics of the JCM. In the same respect it has been shown 
that the measured probability of atomic inversion for specific interaction times turns out to be the 
symmetrically ordered characteristic function [6]. This scheme is closely related to the nonlinear 
atomic homodyne detection [7] in which an atom is coupled to two modes of the field, one acting 



as the signal mode and the other as the local oscillator mode. Furthermore, in 



simple scheme, 
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which can slow down the usual exponential decay of upper state population in an atomic two-level 
system considerably based on an additional intense field with frequency lower than the total decay 
width of the atom, is given. 

Multiphoton single-mode JCM has taken a considerable interest in the literature, e.g. 

For instance, for this model the Heisenberg's equations of motion for the atomic energy 
rave been exactly sloved Its phase variance can exhibit RCP about the long-time 



operator 
behavior 



Moreover, the investigation of this model against squeezed light has shown that the 



111- 



atomic inversion can display RCP for general squeezed input but not for squeezed vacuum 
The analysis of the the model against superposition of squeezed displaced number states, i.e. the 
most general case, is given in 



12j. The multimode version of the JCM has been investigated, in 
particular, the two-mode JCM, e.g. . The most important result related to the two-mode 

JCM is that the atomic inversion exhibits revival-collapse pattern as well as secondary revivals, 
which are independent of the intensities of the initial modes [14]. Moreover, we can mention [lj| 
in which the hamiltonian for the multiphoton multimode JCM has been derived from the first- 
principle. Nevertheless, the generalizations of the JCM as a nonlinear version in both bosonic and 
fermonic variables is given in [16l |. where the exact wavefunction and energy levels are calculated. 

Quasiprobability distribution functions are very useful tools in quantum mechanics since they 
can be used in the calculation of the correlation function of operators as classical-like integrals 
and in the transition to the classical physics. There are three types of such functions, namely, the 
Wigner W, the Husimi Q and the Glauber P functions [17]. These functions are not real probability 
function owing to the position-momentum uncertainty principle. Actually, the W function plays an 
exceptional role among all quasiprobability distributions for several reasons: It contains complete 
information about the state of the system (, i.e. it carries the same information as the density 
operator or as the corresponding wave function). It provides proper marginal distributions for 
individual phase-space variables. It can be used to evaluate the symmetrically-ordered moments 
for the operators of the system. It is sensitive to the interference in phase space and consequently 
it provides a clear prediction to the possible occurrence of the nonclassical effects of the quantum 
mechanical system. In this respect, the W function can be used to analyse the decoherence of 
the quantum system, i.e. the process that limits the appearance of quantum effects and turns 
them into classical phenomena Q, [ljj. It is worth mentioning that the decoherence is useful 



for applications which 
quantum computation 



require keeping coherence in mesoscopic or macroscopic systems such as 



201 ] . Finally, the W function can be determined from the knowledge of the 



complete set of moments of system operators 
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For the single- mode J CM with field prepared initially in coherent light it has been shown that 



there is a relation 
RCP in (a z (T)) 



aetween the behavior of the Q distribution function and the occurrence of the 
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251 ] . For instance, the collapse of the Rabi oscillations in the evolution of 
(<7 2 (T)) is reflected in the behavior of the Q function as the splitting of the initial shifted Gaussian 
distribution into two distributions, which counter-rotate on a circle in the complex plane of the 
distribution. However, the revivals in (a z (T)) correspond to the collision of the two peaks of the Q 
function to produce a single-peak distribution, which is similar to the initial one. It is worthwhile 
mentioning that such relation between the Q function and (a z (T)) for JCM is remarkable only 
when the amplitude of the initial coherent light is very large. Additionally, the comparison between 
(<r 2 (T)) and the Q function has to be performed at the same specific values of the interaction time. 
In this paper we give a new relation, which shows that the information stored in (a z (T)) can be 
obtained from the evolution of the W function at the phase space origin (WOP). This relation 
depends on both the type of the initial state of the optical cavity field, the values of the transition 
parameters and the number of modes interacting with the two-level atom. The motivation of the 
work is two-fold: 

(i) (a z (T)) can be measured using techniques similar to those used for the W function. 

(ii) (a z (T)) can be used to provide information on the nonclassicality of the bosonic system. 
Actually, these are novel results. 

In this paper we consider the interaction of multiphoton N modes of the electromagnetic field 
with a two-level atom in terms of the multimode multiphoton Jaynes-Cummings model JCM. The 
hamiltonian controlling the system is given in the framework of rotating wave approximation. We 
also consider the optical cavity modes are initially prepared in general pure quantum states. For 
this system we seek the relation between the evolution of (a z (T)) and WOP. This will be done in the 
following organization: In section 2 we give the basic relations and equations used throughout the 
paper. In section 3 we discuss the main results as well as we shed the light on how one can measure 
(<7 Z (T)) using techniques similar to those used for the W function. Conclusions and remarks are 
summarized in section 4. 



II. BASIC RELATIONS AND EQUATIONS 

In this section we give the basic relations and equations, which enable us to justify the relation- 
ship between (a z (T)) and the W function for JCM. Firstly, within the dipole and rotating wave 
approximation (RWA) the general form of an idealized hamiltonian, which describes the interaction 
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of multiphoton N modes cavity field with a two-level atom ( JCM) is [la . |26[ ] 

/j- N N N 

— = ^ + ^aVz + A(<T + af + a- a] fcj ), (1) 



where the jth mode is designated by aj (aj) the usual photon annihilation (creation) operator, 
the frequency ojj and the transition parameter kj . &± and a z are the Pauli spin operators describing 
the atomic system, u a is the atomic transition frequency and A is the atom-field coupling constant. 

The hamiltonian can be written as the sum of the two operators: 

N 



3=1 

r ,,„. 



A A ( 2 ) 



C 2 = Aa z + A(«7+ J] af + &- U «f 

3=1 3=1 

where 

N 

£i = fejWj, A = u a - ei (3) 

3=1 

and A is the detuning parameter. Based on the standard commutation rules for the bosonic and 
Pauli operators, it is easy to prove that C\ and C 2 are constants of motion and also they commute 
with each other. In the interaction picture the unitary evolution operator takes the form 

U I (T,0)=exp(-qC 2 ) 

~~ 2^ n! 

n=0 

_ f (-^c 2 ) 2 " , f (zjF^ 

- 2^ (2n)! 2_^ (2n+l)! W 
n=0 ; n=0 V ; 

_ ^ (-1)"(TP) 2 " _ J_ ^ (-l) n (Tv) an + 1 ^ 

— V2n}\ \j> Lj (2n+l)! U2 
ra=0 n=0 



COS(T^) - I \ ' C 2 , 



where 



N N 



T = At, v 2 = (j) 2 + a-o+ H af j af + a+d- \[ a/ af j . (5) 

3=1 3=1 

For making the analysis quite general, we assume that the jth mode is initially prepared in a 
general pure quantum state given by 

oo 

HM) = E ( 6 ) 
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(i) 

where Cn/ represents the probability amplitude for the state under consideration such that 

oo , 

\Cn/\ 2 = 1- We suppose that the atom is initially prepared in the excited atomic state 
|+). Therefore, the initial state of the atom-field system can be expressed as 

|*(0)) = |Vi(0)) (g) l^j(O)) ® • • • <g> I^at(O)} <g> |+) 

(7) 

oo 

= F(ni,n 2 ,--- ,n N )\+,m,n 2 , ■ ■ ■ ,n N ), 
n=0 

where the vector notation in the index means that we have N summations, i.e. n = 
{n\, n 2 , • • • , Un), and the distribution F(m,n 2 , ■ ■ ■ , njv) reads 

N 

F(m,n 2 ,--- ,n N ) = l[C$. (8) 
i=i 

From (jlj) and ([7]) one can easily obtain the dynamical wave function for the system in the 
interaction picture as 

|*(T)) = tfj(r,o)|tt(o)> 



X] F(n 1 ,n 2 ,--- ,n N ) Gi(n 1} n 2 , ■ ■ ■ ,n N ,T)\+,m,n 2 , ■ ■ ■ ,n N ) (g) 

n=0 L 



where 



-iG 2 (nx,n 2 , ■ ■ ■ ,n N ,T)\-,m + h,--- ,n N + k N ) 



h(ni,n 2 ,--- ,n N ;ki,k 2 ,- ■ ■ ,k N )= f] ("^) ! + (A^ 



Gi(ni,n 2 , • • • ,n N ,T) = cos (r^h(n 1 ,n 2 , - • • ,n N ; k t , k 2 , ■ ■ ■ ,k N )^j 



_ ^ sin ^T-^/h(ni,n 2 ,--- ,n N ;ki,k 2 ,--- ,k N )J 
A y/h(ni,n2,— ,n N ;ki,k 2 ,— ,k N ) 



(10) 



sin(T^/h(ni,n 2 ,- .n^fci.fca,- ,feiv)) (n,+fc 7 -)! 

G 2 ni,n 2 ,--- ,n7v,T = y J H J ■ 

yh(ni,n2,— ,n N ;ki,k 2 ,~- ,k N ) \ j=l J ' 

The atomic inversion associated with Q is 

oo 

(<7 Z (T)) = ^|F(m ! n 2j --- ,nAr)| 2 [|G 1 (rii,ri 2 ,- ■ ■ , n N , T)\ 2 - \G 2 (m, n 2 , ■ ■ ■ ,n N ,T)\ 2 }. (11) 

n=0 
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For reasons that will be clear shortly we write down the different forms for the W function. 
The basis of the W function for any quantum mechanical system is the W function of the number 
state | re) having the form 

( — I)" „2 „2\ T /o„2 i o 2\ 



W n (q,p) 



7T 



■exp(-g 2 -^)L n (2^ + 2^), 



(12) 



where L n (.) is the Laguerre polynomial of order re. Also the marginal position probability distri- 
bution for the number state Ire) can be obtained from (1121) as 



(13) 



P(q) = J W n (q,p)dp 



~~ v^F 2"n! 2 >■> 

where H n (.) is the Hermite polynomial of order n. The corresponding form of the marginal mo- 



mentum probability distribution is the same as (|13p but q should be replaced by p. The A^-mode 
dynamical W function can be defined up to a constant prefactor as [itJ 



W(J3,T) = TV 



N 



(14) 



where /3 = (A, /%,••• ,#v) = (<7i, <72, • • • , QN)Pl,P2, • • • ,Piv) since /?,• = qj + ipj- p(T) is the density 
matrix for the system under consideration and -D(/?) is the multimode displacement operator having 
the form 

"~ AT 



£>(/?) = exp 



(15) 



At the phase space origin (, i.e. = 0) the formula (fill) reduces to 



PF(0,r) = Tr p^exp UTT^atojj . (16) 

Formula (|16p indicates that the main contribution for the W function at the phase space ori- 
gin is resulted from the diagonal part of the density matrix of the quantum mechanical system. 
Comparing this situation with that of the atomic inversion one can conclude that there is a clear 
relationship between the evolution of WOP and the corresponding atomic inversion. 
Now on substituting ([9]) into (|16p and carrying out the expectation value we arrive at 



W(0,T)= £ \F( ni ,n 2 ,--- ,n^)| 2 (-l) ni+ * a+ - +nw {|G ! i(ni,n a ,-- - ,n N ,T)\ 2 

n=0 1 



+ (_l)fci+-+fcw|Gf 2 ( 



rei,n 2 , • • • ,n N , 



T)\ 2 }. 



(17) 
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Expression (|17p and its consequences are the main results of the paper. Specifically, we show 
that for particular types of the initial states and particular values of the transition parameters kj, 
expression (fT7|) coincides with that of the corresponding atomic inversion. This will be discussed 
in the following section. 

We proceed by connecting the present results with those of the homodyne tomography tech- 
nique. Therefore, we give the mathematical relation between the iV-mode W function and the 
corresponding distribution function pr(q\,--- ,qjy,0i,--- ,0n,T) (, i.e. Radon transformation). 
Such relation is just the generalization of the single-mode case and can be expressed as 

oo oo 

pr{qi,- ■ ■ ,qN,0ir ■ ■ ,&N,T) = f dpi--- f dp N 

— oo —oo 

(18) 

xW(q\ cos #i — pi sin 0±, qi sin 0\ + p\ cos 6\, ■ ■ ■ ,qjy cos On — pjy sin On, qN sin On + Pn cos 9n,T). 

In (|18p we have assumed that N modes can be delivered to N separate ideal balance homodyne 
detectors. At the phase space origin, i.e. qj = 0,0j = 0,j = 1, •••iV, the formula (fT8j) becomes 
phase independent and reduces to 

oo oo 

pr(0,T)= J dpx--- f dp N W(0, Pl ,0,p2,--- ,0,PN,T), (19) 

— oo — oo 

where W(0,px,Q,p2, ■ ■ ■ ,0,pn,T) is the diagonal part of the W function, which is phase indepen- 
dent. On using (fT2j) and (fT3|) one can easily deduce pr(0, T)(= P(0,T)) given by (fT9j) for the state 
vector Q as 



oo N h 2 (0) 

P(0,T) = Z \F( ni ,n 2 ,--- ^^iMlGxK,^,--. ,n N ,T)\ 2 ft 

n=0 L j=l 3 ' 



(20) 



N H 2 , , (0) 



3=1 



III. MAIN RESULTS 



In this section we discuss two issues: (i) We investigate the results by making a comparative 
study among the behavior of (a z (T)),W(0,T) and P(0, T), for the system under consideration, 
(ii) We argue how one can measure the atomic inversion via techniques similar to those used for 
the W function. 
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A. Investigation of the results 

As we mentioned above we investigate the behavior of the quantities W(0, T) and P(0, T), and 
then compare such behavior with that of the corresponding (a z (T)). 

We start the discussion with W(0, T), which is given by (|17p . It is obvious that when m + n<i + 
■ • • + Wjv is even and k\ + ki + • • • + kjq is odd W(0, T) is identical with the atomic inversion of 
the system (c.f. (Ilip ). In other words, the atomic inversion can be used to provide information on 
the nonclassicality of the dynamical bosonic system. For instance, when the evolution of (a z (T)) 
displays negative values, the JCM can exhibit nonclassical effects. Nevertheless, when k\ + ki + 
• • • + kx is even number regardless of the value of n\ + ri2 + • • • + n^, the expression (fTT|) becomes 
time independent (, i.e. W(0, T) is localized) and can be factorized in the following sense 

N 

W(0,T) = Y[Wj(p,0), (21) 
i=i 

where Wj (0,0) is the initial value of the W function of the jth mode at the phase space origin 
having the form 

oo 

W 3 (0,0)= £(-lHC#| 2 . (22) 

n j= 

The expression f)21 1) can be obtained, e.g., when the number of modes iV is even and the transition 
parameters are symmetric, i.e. k± = k% = ■ ■ ■ = k]y. Moreover, expression (|2ip indicates that if the 
initial W function of only one of the modes has a negative value at the phase space origin whereas 
those of the others are positive, the system can provide nonclassical effects. This is a sufficient but 
not necessary condition. 

Now we give a closer look at the behavior of the W(0,T) for the single-mode case, i.e. N = 1. 
In this case expression (fT7|) reduces to 

oo 

W(0,T)= IC^W-ir^G^n 2 + (-l)*i|G 2 („ 1) T)| 2 }. (23) 

ni=0 

For odd transition parameter and initial even (odd) parity states, e.g. even (odd) coherent states, 
(|23|) gives 

W±(0,T) = ±(a z (T)), (24) 

where "+" and "-" signs denote even and odd parity states, respectively. For A = expression 
(|24j) indicates that when the initial intensity of the radiation field is weak the system can exhibit 
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FIG. 1: The evolution of the atomic inversion (a z (T)) (a) and the W function W(0,T) (b) against the 
scaled time T for the single-mode case with k\ — 1 and when the field and atom are initially prepared in 
the coherent state (with \a\ = 8) and atomic excited state, respectively. 

nonclassical effects periodically. Nevertheless, in the strong-intensity regime one has W(0, T) ~ 
0), which is associated with the occurrence of the collapse (revival) in (<t 2 (T)). Consequently 
the W function exhibits nonclassical interference at the phase space origin only in the course of 
the revival times, i.e. the nonclassical effects most probable occur in the course of the revival time. 
However for the non-parity states the locations (in the interaction time domain) of collapses and 
revivals occurring in W(0,T) are interchanged compared to those in (a z (T)). This agrees with 
the fact that the JCM generates Schrodinger-cat states in the course of the collapse time [121 ]. We 
proceed by investigating the behavior of the W(0, T) for the standard JCM, i.e. ki = 1, A = and 
the field is initially prepared in coherent light with amplitude \a\. In this case (I23p reduces to 

W(0,T) =exp(-|a| 2 ) £ J£Li(-l)ni cos(2T^hTT) 

(25) 

= exp(-|a| 2 ) £ HLicoBCZTv^r + T + niTr). 
m=0 

Expression (|25p is identical with that of the corresponding atomic inversion but with additional 
factor, which is (— l) ni . This factor is responsible for the interchange of the "locations" of collapses 
and revivals occurring in W(0,T) compared to those exhibited in (a z (T)), as we mentioned above. 
This is remarkable in Figs. 1(a) and (b) where we have plotted (a z (T)) and W(0,T), respectively, 
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FIG. 2: The evolution of the W(0, r) against the shifted-scaled time r for the same situation as in Fig. 
1(b). 



for given values of the interaction parameters. 

The evolution of the atomic inversion (a z (T)) (a) and thethe W function W(0,T) (b) against 
the scaled time T for the single-mode case with k\ = 1 and when the field and atom are initially 
prepared in the coherent state (with \a\ =8) and atomic excited state, respectively. 

Also this can be emphasized by deducing the asymptotic form for f)25 j) in the strong-intensity 
regime (, i.e. \a\ is large). By means of the harmonic approximation technique (see equation 
(la) in the appendix) and after straightforward calculations (|25p can be expressed as 



W(0, T) = exp 



-2n cos 2 (— ) 
2n 



cos 



1 T 
Tin + — ) — nsin( — ) 
n n 



(26) 



where fi = \a\. Expression (I26p is similar to that of {a z (T)) except cos(.) in the exponent should 
be replaced by sin(.). As a result of this fact the envelope function in (|26p gives its maximum 
value at T = 7rn, whereas, that of (a z (T)) is maximum at T = 2nn. From these arguments and 
information displayed in Figs. 1 one can conclude that W(0, T) can give similar information on the 
corresponding atomic inversion provided that the interaction time T is replaced by r = T + irn. 
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FIG. 3: The probability distribution measured using homodyne tomography corresponding to the evolution 
of W(0, T) against both the scaled time T (a) and shifted-scaled time r (b) for the same situation as figure 
1(b) and figure 2, respectively. 

The behavior associated with this situation is given in Fig. 2. Comparison between Fig. 1(a) and 
Fig. 2 is instructive. 

Now we turn the attention to the formula (|20p . which is related to the homodyne tomography. 
Firstly, it is worth reminding that the properties of the Hermite polynomial provide H2 n +i (0) = 
and H2 n (0) ^ 0. These facts make the argument related to P(0, T) different from that given for 
W(0, T). For instance, when one of the modes is initially prepared in odd parity states (, e.g. odd 
coherent states) and the associated transition parameter with this mode is an even number then 
-f (0, T) = 0, however, this is not the case of the corresponding W function, where W(0, T) ^ 0. 
Also one can easily recognized that P(0, T) ^ for different cases, e.g., when all modes are initially 
in even parity (non-parity) states regardless of the values of kj or when all modes are initially in 
odd parity states provided that the transition parameters are odd numbers. 

Similar to the treatment given for the W function we investigate P(0, T) of the single-mode 
case when the field is initially in coherent state, k\ = 1 and A = 0. Therefore, relation (|20p can 
be expressed as 
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Laser source 




FIG. 4: The proposed setup for detecting (a z (T)) of the single-mode JCM when the field is initially prepared 
in the coherent light pumped by the laser source. BS and PD denote beam splitter and photodetector, 
respectively. The two-level excited atom is localized in the microcavity as indicated. The initial state, 
detected state and local-oscillator state are denoted by \ip(Q)), \4>{T)) and \oll), respectively. 



P(0,T) = ± £ P(m; 



711=0 



+ 



H^(0) 
2"lni! 



H£ +i(0) 
2"i+ 1 (ni+l)! 



gyp) h^ i+1 (o) 

2"mi! + 2"i +1 (ni+l)! 



cos(2r^TTT)}, 



(27) 



where P(ni) is the photon number distribution for coherent states. In the strong-intensity regime 
the asymptotic form for (|27p. which is corresponding to (|26p. is 



P(0,T) = iexp(-|a| 2 ){l (|a| 2 ) + I 1 (|a| 2 )+ 



(28) 



Re{exp[iT(n + I)] [l (|a| 2 exp(i|)) - Ii(|c| 2 exp(ig))] }}, 

where Io(-) and Ii(.) are the modified Bessel functions of the first kind of order zero and one, 
respectively. The derivation for (128p is given in the appendix. Expression (128D is periodic with 
period 2-7rn where n is an integer. Additionally, P(0, T) gives its maximum values around T = mn-zr 
whenever m is an odd integer. We have plotted (|27j) in figures (3) against both the scaled time 
T (a) and the shifted-scaled time r (b) for the same situations as those for figure 1(b) and figure 
2, respectively. Comparison between Fig. 1(b) and Fig. 3(a) as well as Fig. 2 and Fig. 3(b) 
shows that RCP in (a z (T)) can be observed via homodyne tomography. Of course, the evolution 
of (a z (T)} and P(0, T) possess different scales. 

This problem can be solved easily by using (|28p and the information displayed in Fig. 3(b). For 
instance, in the strong-intensity regime we can adopt the following relation 



Pn 



[P(0,r)-P(0,0)], 



(29) 
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where 

P(0,0) =exp(-|a| 2 )I (|a| 2 ), 

(30) 

Pmax = [I (|a| 2 ) + Ii(|a| 2 )]exp(-|a| 2 ), 

where the subscript max stands for the maximum value of P(0,T). The explicit form for P ma x 
can be obtained by analysing the behavior of P(0, T) around T = nn. Finally, the origin of the 
prefactor (, i.e. 1/P m ax) m f|29|) can be understood as follows. P(0,T) carries complete information 
on (a z (T)) when it has not only similar behavior but also the same amplitude as (a z (T)). Thus 
we seek an amplification factor, say, fi (, i.e. /iP(0,T) = (a z (T))) such that 

I^Pmax = | (<5"z(P))max| = 1- (31) 

B. Observation and measurement 

In the first part of this section we have shown that under certain conditions there is a direct 
relation between (a z (T)) and W(0,T) for JCM. Such relation indicates that the atomic inversion 
can be detected via techniques similar to those used for the W function. This will be discussed in 
the following, in particular, for the single-mode JCM. 

As it is well known there are different schemes proposed for measuring the W function, which 



are: Photon counting experiment 



QED and ion trap s 



measurements 



31 



29 



28 1, using simple experiment similar to that used in cavity 
30j, and tomographic reconstruction from data obtained in homodyne 
321 ] . Here we argue how these techniques can be used for measuring (a z (T)). 
(i) Measurement of the (a z (T)) using photon counting experiment: Photon counting method is 
based on the fact that the single-mode W function at the origin of the phase space can directly 



281 ] . Hence (a z (T)) can be detected via this 



be measured by a photodetector facing this mode 
technique in the following sense (see Fig. 4) . Generally, a mode prepared in a coherent state, which 
is pumped by a laser source, interacts firstly with the two-level excited localized atom — localized 



and/or very slow atom can be prepared by means of, e.g., laser-cooling technique 33j — hence the 
outgoing field is superimposed by a strong local-oscillator mode (|«l)) via beam splitter whose 
transmissivity (reflectivity) is high (low). The measurement has to be performed only on one of 
the output ports of the beam splitter via the photodetector. It is worthwhile mentioning that the 
one-port measurement for the beam splitter can be performed u sing a conditional measurement 
technique in which no photons are measured in the free port, e.g., [34]. We proceed by considering 
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a perfect photodetector, i.e. its efficiency is unity. The probability P{n\,T) of the registration of 
n\ photons at time T in the detector is given by 

P(n 1 ,T) = (:^-exp(-A):) J (32) 

where :: stands for normally-ordered operator, angle brackets mean expectation value (which is 
calculated in the framework of Schrodinger picture) and A is the operator of the integrated flux of 
light onto the surface of the detector having the form 

A = {Vfa ] - VI - fa* L )(Vfa - Vl - fa L ), (33) 

where f is the transmissivity power of the beam splitter. The count statistics determined in the 
experiment is used to compute the following count generating function 

oo 

G(a,T)= Z (-l) ni P(ni,r) 

(34) 

= (: exp(-2A) :). 

The second line in (I34|) is obtained by substituting (I32p in the first line of this equation. It is 
evident that for A = 0, W(0,T) given by ([23]) can be expressed in the form flMJ) , where, in this 
case, we have 

P(m,T) = \C^\ 2 [cos 2 (T^hinuh)) + sin 2 (T ^ h(n u -h))j . (35) 

We proceed from the definition of the W function (I16p and (134)) the W function for the detected 
mode is proportional to s = 1 — 1/f ordered quasidistribution function of the mode entering the 
beam splitter (which is not included in Fig. 4) as 



G(a,T) = lw ^^«L,^,Tj . (36) 

When the transmissivity of the beam splitter is near one the scanned quasidistribution is W(0, T), 
i.e. the atomic inversion is detected. Here we assume that the interaction time in the microcavity 
and the detection time in the photodetector are equal since we are interested only in showing that 
(<t z (T)) can be measured via photon counting technique. Nevertheless, in the realistic situation 
there is a time delay between the interaction and the detection processes. This problem can be 
solved by obtaining information on both the detector-microcavity distance and the velocity of the 
radiation field. 



15 



(ii) Measurement of the (<7 2 (T)) using cavity QED: In this technique, probes of the two-level 
atoms interact dispersively with the field under consideration-the source of the field in this case 
is the JCM to which we would like to measure (<r 2 (T) )-causing a phase shift to the atomic wave 
function, which is proportional to the photon number. This phase shift can be revealed by Ramsey 
atomic interferometer 35J. Under certain conditions the final state of the atom measures the field 



parity at the phase space origin, i.e. W(0, T). The experimental setup related to this proposal is 



given in [29( (see Fig. 1 in 



291 ]) but with additional arrangements. For instance, in this setup we 



have to block the microwave generator, which makes displacement for the field under consideration 
since the interest is focused only on the behavior of the W function at the phase space origin. 

(iii) Measurement of the (a z (T)) using homodyne tomography arrangement: Generally single- 
mode homodyne tomography is based on the set of distributions pr(qi,6i) measured by homodyne 
detection, i.e. the field to be measured beats with the local oscillator in a homodyne arrange- 
ment [3^ . [37}]. Once obtained pr(qi,0i) the W function can be reconstructed via inverse Radon 
transformation. At the phase space origin this relation reads 



OO 



W(0,T) = ±- J dC J dr/MpKCTjexp^C), (37) 

—00 —00 

where pr(£,T) = (£\p(T)\Q. The value of the relative phase between the local oscillator and the 
signal field, which is assumed to be the field outgoing from the JCM microcavity, is zero. This can 
be arranged by moving a mirror on a piezoelectric translator [38]. It is worthwhile mentioning that 
pr((, T) has been measured via this technique, e.g., in Q. Also investigation for the random-phase 
states using homodyne tomography technique is given in [39(] . 

Finally we conclude that the techniques (i) and (ii) can lead to a direct measurement of {a z {T)). 
Also, they do not involve inversion algorithm and hence they should be much less sensitive to 
experimental errors than the tomographical technique. Furthermore, they can be, in principle, 
applied to the atomic inversion of the multimode JCM where the attention has to be focused on 
the measurement of the (a z (T)) in its entangled form. 



IV. CONCLUSIONS AND REMARKS 



In this paper we have discussed the relation between both the evolution of the atomic inversion 
and the corresponding W function for multimode multiphoton JCM. We have shown that under 
certain conditions there is a direct relation between these two quantities, which is valid for resonance 
and off-resonance cases. Such relation suggests that the nonclassical effects stored in the radiation 
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field can be noticed through the behavior of {a z (T)). Furthermore, based on this relation we have 
discussed the possibility of detecting (a z (T)) for single-mode case using techniques similar to those 
applied to the W function. 

The results given in this paper are valid to any JCM hamiltonian provided that it has been de- 
duced in the framework of rotating wave approximation. More illustratively, if in (pQ) the operators 

Y[ o,J are replaced by a J 1 a^ 2 a^ 3 a\ 4 • • • , i.e. some of the annihilation operators are replaced by 
3=1 



creation ones, the relation (|17p will not be affected owing to the fact 

(n - k\ exp(iirdt a)\n - k) = (n + k\ exp(iirat a)\n + k) = (-l) n+k . (38) 

Moreover, the results are independent of the type of the initial atomic state (, i.e. if the atom is in 
the excited state, ground state or atomic superposition state). Basically expression (|T7|) depends 
on the Fock state basis of the dynamical wave function, e.g. for the single-mode case it depends on 
\n) ,\n + k) and \n — k). Consequently, such relation still exists when Jaynes-Cummings hamiltonian 



includes Kerr nonlinearity 40(] , Stark effect 41[ , intensity dependence [42j and atomic motion 43J] . 
Furthermore, among all quasiprobability distribution functions such relation, i.e. (|17p . exists only 
for the W function. For instance, the evolution of the Q function at the phase space origin for 
single-mode JCM with initial coherent light and A = is given by 

Q(0, T) = exp(-2|a| 2 ) cos 2 ! (T '\/h\) . (39) 

It is obvious that in the strong-intensity regime Q(0,T) — * 0. 

Finally, as it is well known, the W function is a global quantity which characterizes the full 
quantum state. Additionally, dealing with the W function at an isolated single point causes a 
difficulty in finding a proper normalization. Nevertheless, throughout the paper we have focused 
the attention on the behavior of WOP and its consistence with the behavior of (a z (T)). Therefore 
the normalization has no effect on the dynamical behavior of the system. 



Appendix 

In this appendix we derive the asymptotic form f|28j) for P(0,T). It is worth remembering that 
in the strong-intensity regime (, i.e. \a\ >> 1) the argument of cos(.) in (|27p can be expressed as 
2J: 

v^+T = y/(fi) + n + 1 - (h) ~ ±(n + \ + |), (la) 
where n = \J (n(0)). 
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Now we show how the different summations in (1271) can be evaluated. 



CO . .o 2-7T CO ( _3_\n CO ( cf_\m 

E (J^ H n(0) = h I E 4^(0) £ ^-H m (0)#, (2a) 

n=0 "=0 m=0 

where a = |a| exp(i<^>). By means of the generating function of the Hermite polynomial the sum- 
mations on the right hand side of (2a) can be carried out as 

OO . |2n 2-7T 

E w^ H ^°) = h I exp[-|a| 2 cos(2< ) 6)]# = I (|a| 2 ), (26) 
n=0 v o 

where Io(-) is the modified Bessel function of the first kind of order zero. The second summation 
we would like to evaluate is the following 

\a\ 2n r2 fn s _ ^ |a| a(m-i) 2 
n!(n+l)!2"+l n n+lV U J ~ 2^ ( m ~l)!m!2 m n ml U J> W 
n=0 ( ' m=0 V ; 

where the factorial — l! is oo. Summation (3a) can be reformulated as 

oo i^|2n 2 ^ < ^ |ap m 2 

E n!(n+l)!2«+l H n+l(°) = dR 7 E ( m \)'l 2 ™ H m( )- ( 3& ) 
n=0 m=0 

Using (2b), the right hand side of (3b) gives 

^I (|a| 2 ) = I 1 (|a| 2 ). (3c) 

The terms including cos(.) in the strong-intensity regime can be written as 

E 7^rH 2 (0)cos (T(n + ± + f)) = Re jexp[iT(n + I)] £ t&^H 2 (0) exp(zT|)}. (4a) 

The summation on the right hand side of (4a) can be evaluated using procedures as those given 

above leading that 

E t^H 2 (0)cos (Tin + i + f)) = Re{exp[iT(n + I)]I (|a| 2 exp(i|))}. (46) 

n=0 V ' ; L J 

Similarly the last summation in (|27[) can be performed. 
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